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ABSTRACT 

A plausible model for the Crab Nebula is one in which a particle dominated, highly 
relativistic wind from the pulsar passes through a shock front in which the particles at- 
tain a power law energy distribution. The electrons and positrons lose energy radiating 
synchrotron emission. Here, a one zone version of the model is developed and applied 
to observations of X-ray pulsar nebulae. Efficient conversion of pulsar spin-down power 
to X-ray luminosity is expected if the observed electrons are in the synchrotron cooling 
regime and their energy spectrum is similar to that in the Crab Nebula, which has 
N(j) oc 7~ 2 ' 2 , where 7 is the particle Lorentz factor. In this case, the relation between 
X-ray luminosity and pulsar spin-down power depends only weakly on the model pa- 
rameters. The dependence is stronger for a steeper particle spectrum, as appears to be 
present in N157B, and the efficiency of X-ray production can be lower. If the electrons 
are not in the cooling regime, the X-ray luminosity can be low, as appears to be the 
case for the compact nebula around the Vela pulsar. Slow cooling is likely to be caused 
by a sub-equipartition magnetic field in the radiating region. Observations can place 
constraints on the uncertain physics of relativistic MHD shocks. The model is related 
to those developed for gamma-ray burst afterglows. 

Subject headings: ISM: supernova remnants — shock waves — stars: pulsars: general 
— X-rays: general 



1. INTRODUCTION 

Beginning with the Crab Nebula, extended sources of X-ray synchrotron emission have been 
observed around a number of pulsars. The power for the relativistic particles and magnetic fields is 
thought to be the spin-down power of the pulsar, E. In cases where E can be measured, there is a 
relation between the X-ray luminosity, L x , and E. On the basis of Einstein data, Seward & Wang 
(1998) found a relation of the form logL x = 1.39 log E — 16.6, where L x is the X-ray luminosity in 
the 0.2-4 keV range and E is the pulsar spin-down power in ergs s _1 . The X-ray luminosity of the 
central pulsar is included in addition to the extended nebula. On the basis of ASCA data, Kawai, 
Tamura, & Shibata (1998) found the relation logL^ = (33.42 ± 0.20) + (1.27 ± 0.17) log (£/10 36 ), 
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where L x is now the nebular luminosity in the 1-10 keV range. Expressions like these have been 
widely used to estimate E in X-ray synchrotron nebulae where a pulsar has not been detected. 

The objects considered by Kawai et al. (1998) cover a range of E of 4x 10 33 — 5x 10 38 erg s _1 and 
include symmetric objects like the Crab Nebula, which may be interacting with supernova ejecta, 
and the CTB 80, which appears to be a bow shock nebula interacting with a supernova remnant. 
A compelling model for the optical/X-ray properties of the Crab Nebula was developed by Rees & 
Gunn (1974), Kundt & Krotscheck (1980), and Kennel & Coroniti (1984a, 1984b; hereafter KC84a 
and KC84b). In this model, the central pulsar generates a highly relativistic, particle dominated 
wind that passes through a shock front and decelerates to match the expansion velocity set by the 
outer nebula. The electrons and positrons in the wind acquire a power law energy spectrum in the 
shock front and radiate synchrotron emission in the downstream region. KC84b showed that this 
model is able to account for the basic high energy properties of the Crab, including the position of 
the wisps, the spectrum of the nebula, and the size of the nebula at different wavelengths. They 
were able to assume a steady state flow because the particle lifetimes to synchrotron losses are less 
than the age of the Crab Nebula. The model parameters are the Lorentz factor of the wind, "y w , 
the energy spectral index of the particles accelerated at the shock front, p, and the magnetization 
parameter of the wind, a, which determines the magnetic field strength in the extended nebula. 

Here, I develop a simplified version of the shocked wind model for the Crab Nebula which 
captures the essential features for estimating the X-ray luminosity and spectrum (§ 2). The model 
can be compared to the properties of other pulsar nebulae in order to investigate the implications 
for the model parameters (§ 3). 



2. X-RAY EMISSION MODELS 

The model developed here is similar to that of KC84a and KC84b, but is a one zone model that 
does not attempt to account for the spatial structure of the nebula. The pulsar spin-down energy 
is assumed to go into a highly relativistic wind with Lorentz factor 7^. As noted by KC84a, the 
magnetization parameter in the wind, a, must be small (< 0.1) in order to have a high efficiency 
of wind power into synchrotron radiation. With these approximations, we have 

E rj 47mi7^r;W 3 , (1) 

where n\ is the proper density in the wind just before the shock front, r s is the shock wave radius, 
m is the electron mass, and c is the speed of light. 

At the shock front, the wind particles acquire a power law energy distribution of the form 
iV(7) oc 7~ p for 7 > 7 m , where n = J ^(7)^7 is the density of relativistic particles and 7 is the 
particle Lorentz factor. I assume that p > 2 and that the particle spectrum extends to higher 
energies than those of the X-ray emitting particles. In the first model, the particles responsible 
for X-ray emission rapidly lose their energy to synchrotron radiation so there is a steady state. 
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The balance between injection from the shock front and synchrotron losses leads to a constant 
number of emitting particles. The basic simplification made here is that the emitting region can 
be treated as one zone. The energy density in the zone is approximately determined by the shock 
jump conditions. The particle pressure behind a highly relativistic shock wave is 2nimc 2 7 2 J /3 
and the dynamic pressure is nirac 2 7^,/3 (KC84a). The flow decelerates downstream from the 
shock front, so that nimc 2 ^ is an estimate of the pressure and 3nimc 2 ^ is an estimate of the 
energy density. I assume that the energy density in the emitting region is divided between a 
fraction e e in particles and a fraction e_e in the magnetic field. These parameters, which sum to 
unity, replace the wind magnetization parameter, a, of KC84a. In the flow model of KC84a, a is 
small but the compressive flow downstream from the shock front builds up the magnetic field to 
close to equipartition. The behavior of the magnetic field in the downstream region is not well 
understood. The toroidal field structure of the KC84a model is likely to be unstable (Begelman 
1998), and polarization observations of the synchrotron emission in the Crab Nebula show that the 
field structure is complex. In the relativistic blast wave model for gamma-ray burst afterglows, some 
form of magnetic field enhancement in the postshock region is required and it is often characterized 
by an efficiency factor €b like that defined above (e.g., Sari, Piran, & Narayan 1998). 



With these assumptions, the magnetic field in the emitting region is 

1/2 

^ r 2 s c ) 



The postshock density in the shock frame is 3j w ni for a strong relativistic shock, so that the 
postshock particle distribution function is iV( 7 ) = (jp — l)3 7u ,7m _1 ni7 _p . The expression for the 
energy density then leads to 

The number of radiating particles at a particular 7, A/"(7), is determined by a balance between 
the rate at which particles are injected at the shock front, A/" (7), and synchrotron losses. The 
radiation power of an electron with Lorentz factor 7 S> 1 is 

P( 7 ) = ia T c7 2 f^/3i?V, (4) 

where ot is the Thompson cross section and j3 = 1.06 x 1CP 15 cm 3 s _1 . The total rate at which 
particles are injected is E/^wmc 2 ), so that 

AT( 7 ) = (p _ \) 1 P n ;\ lw mc 2 )- l E 1 - p (5) 

for the spectrum determined above. The balance between injection and synchrotron losses can be 
expressed as 

1 dM/3B 2 j 2 



mc 2 $7 



= AA( 7 ) (6) 
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with the solution 

N( 1 )=Y m \l w PB 2 )- 1 E 1 -^ +l \ (7) 
I assume that the electron energy is radiated at its critical frequency 

where q e is the electron charge and £ = 2.80 x 10 6 in cgs units. The luminosity of radiating particles 
in the range 7 to 7 + is f3B 2/ y 2 N('~v), which leads to a luminosity per unit frequency 

1 / o\P~ 1 /ct2\(P — 2 )/4 

L „ = I f Ez| j (*L\ ^WI^2 r7 ip-W&P+W v -vl*. (9) 

This expression is valid provided the X-ray frequency of observation, vx, is greater than the cooling 
frequency v c and u m , corresponding to 7 m , and is less than u max , corresponding to the maximum 
7 in the power law spectrum. The cooling frequency corresponds to the 7 at which the electrons 
are able to radiate their energy in the age of the wind nebula, t, and is 

In the case of a bow shock nebula, the flow time past the nebula may be the relevant value of t 
rather than the total remnant age. The other critical frequency is 

P -2\ 2 e 2 el Uf6e B E\ 1/2 (n) 



p-lj r s \ c J 

The steady state assumption requires that the observing frequency, u D , be > u c . If v c < v < 
v rn , then jv (7) = in eq. (6), which has the solution ^(7) oc 7~ 2 . This yields L v oc u~ l l 2 (see 
also Sari et al. 1998). The constant of proportionality is determined by joining the solution onto 
the v Q > u m solution. 

Observations of the compact nebula around the Vela pulsar indicate v c > vx and a different 
model is needed (see § 3). The Vela case is apparently a bow shock nebula created by the proper 
motion of the pulsar (Markwardt & Ogelman 1998). The nebula age is then t w r s /v p , where 
v p is the transverse pulsar velocity. Because t is less than the time since the birth of the pulsar, 
the assumption of constant E is adequate. The number of radiating particles is ^(7) = Jv(j)t = 
(p — l)7m~ (7w"ic 2 )~ 1 Er s Vp 1 ~f~ p and the luminosity is 

K = 1 (P - XT 1 6^ e P-i € (r+W -(p-i)/2 -1 p^^y^-ip-m (12) 



which yields a lower luminosity than eq. (9) because of the limited number of radiating electrons. 
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3. THE CRAB AND OTHER PULSAR NEBULAE 

Detailed models have been produced for the Crab Nebula, so it provides a point of reference 
for the approximate model developed here. The X-ray spectrum of the Crab is a power law with 
photon index 2.1 (e.g., Pravdo & Serlemitsos 1981), which corresponds to p = 2.2 in the present 
cooling model. With this value of p, equation (9) becomes 

L u = O.OS^e^Vj ^ ^- 11 , (13) 

where cgs units are used. The values r s = 3 x 10 17 cm and B = 5 x 10 38 erg s _1 are determined 
by observations of the Crab Nebula and its pulsar (see KC84b). The efficiency factors, constrained 
by e e + 6b = 1, are not well determined. In the model of KC84a, the particles and magnetic field 
come into approximate equipartition (e e = e# = 0.5), although the field is initially weaker. The 
need for a particle dominated shock wave suggests that e e > 0.5. The quantity e\ 2 e B m is 0.42 in 
equipartition, rises to 0.81 for €b = 0.05, and slowly drops to 0.63 for €b = 10~ 4 ; there is little 
sensitivity to this parameter. I set it to 0.5. The value of ^ w depends on features in the nebular 
spectrum and cannot be determined just from X-ray observations. KC84b use j w = 3 x 10 6 in their 
Fig. 13. When these parameter values are substituted into eq. (13), the result at v = 10 18 Hz is 
vL u = 1.0 x 10 37 erg s _1 , which compares well to the detailed model presented by KC84b in their 
Fig. 13. The model 0.2 - 4 keV luminosity is 2.6 x 10 37 erg s _1 . 

Although the luminosities agree, there are differences between the models. In the present 
one zone model, the synchrotron burn-off leads to a steepening of the frequency spectral index 
by 0.5, the standard result. In the model of KC84b, the particles initially move out through a 
region in which B oc r if the wind magnetization parameter is small. If the synchrotron burn-off 
occurs in this region, the index steepens by (p + 7)/18, which is > 0.5 for p > 2. The additional 
steepening occurs because at a high frequency, the particles typically radiate in a lower magnetic 
field and higher 7 particles are observed. The fact that there are fewer of these particles steepens 
the spectrum. However, for p = 2.2, the steepening is by 0.51, very close to the one zone model. 
For the above Crab Nebula parameters, eq. (10) gives v c = 4.2 x 10 12 Hz for e B = 0.5 and eq. (11) 
gives v m = 1.3 x 10 14 Hz. The fact that u c < u m implies that a spectrum of the form L u oc i/~ - 5 
should occur between these frequencies. Fig. (13) of KC84b shows that such a region does appear 
to be present at optical wavelengths. 

Table 1 gives a summary of well-established X-ray nebulae around pulsars. Kawai et al. (1998) 
list 6 additional nebulae, but Becker et al. (1999) and Pivovaroff, Kaspi, & Gotthelf (2000) do not 
confirm the diffuse emission seen with ASCA in 5 of the cases. The value of E is determined 
from E = where / = 10 45 g cm 2 is the assumed neutron star moment of inertia and is 

the pulsar spin rate. The X-ray luminosities are mostly from the Einstein observatory and cover 
the energies 0.2-4 keV. Several of the photon indices are from ASCA because of its relatively 
broad energy coverage. Two spectral indices are listed for Gil. 2-0. 3; this remnant shows thermal 
and nonthermal emission and the first one fits the lines with Gaussians, while the second uses a 
non-equilibrium ionization model (Vasisht et al. 1996). 
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N158A has properties that are very close to the Crab Nebula, but N157B has a lower value 
of L x /E. The observed photon index for N157B implies p ~ 3 if it is in the cooling regime. With 
p = 3, e e = €b = 0.5, and j w and r s as for the Crab Nebula, eq. (9) yields L x (0.2 — 4 keV) = 
3 x 10 36 erg s _1 . The lower efficiency is naturally produced; the discrepancy with the observed 
luminosity is not a problem because the higher value of p leads to a greater sensitivity to the model 
parameters. 

The compact nebula around the Vela pulsar has an especially low value of L x /E and a relatively 
flat spectrum. The value of L x I have listed here is lower that often quoted because I have taken a 
distance of 250 pc, as advocated by Cha, Sembach, & Danks (1999) and others; even at 500 pc, it is 
a low luminosity object. The properties of the nebula indicate that it is not in the cooling regime, 
i.e. v c > v x for Vela, de Jager, Harding, & Strickman (1996) note that the nebular luminosity 
out to 0.4 MeV is 0.1% E, implying a high value of v c (> 10 20 Hz). The value of r s for Vela 
can be estimated from the size of the compact nebula, which appears to be a bow shock nebula 
(Markwardt & Ogelman 1998). At 250 pc, the radius is 2 x 10 17 cm. The proper motion of the 
pulsar (Ogelman, Koch-Miramond, & Auliere 1989) gives a transverse velocity of 45 km s _1 so 
that t = r s /v p = 4.4 x 10 10 s. The age of the Vela remnant is ~ 10 4 yr, so eq. (12) should be 
applicable and yields L x (0.2 — 4 keV) = 2 x 10 34 (es/10~ ) 0,83 erg s _1 . The non-cooling behavior 
requires v c > 10 20 Hz or 6b < 10 4 . In a similar argument, de Jager et al. (1996) set t = ry/3/c 
as the escape time and thus found a limit on the magnetic field close to equipartition. The model 
here assumes an MHD flow through a bow shock region, which gives a longer residence time for the 
particles. The present model implies considerable variation in es among nebulae; such variation 
has little effect on L v in the cooling case (see eq. [13]). 

CTB 80 appears to be in the same class as the Vela nebula; it has a bow shock structure (Safi- 
Harb, Ogelman, & Finley 1995) and has a relatively hard spectrum and low efficiency of conversion 
to X-ray emission (Table 1). W44 also has a bow shock nature (Harrus, Hughes, & Helfand 1996), 
but the spectrum and X-ray luminosity are too uncertain to draw clear conclusions. 

The other 2 objects in Table 1 have L x /E ~ 0.01. If the nebulae have p ~ 2.2 and rapid cooling 
in the X-ray regime, the weak sensitivity to the parameters in eq. (13) makes it difficult to reduce 
L x /E by the required factor of 5 compared to the Crab. However, MSH 15-52 appears to have a 
flatter X-ray spectrum than the Crab, du Plessis et al. (1995) find that there is steepening of the 
X-ray spectrum, consistent with p = 2.2 and a synchrotron cooling break at ~ 6 keV. A cooling 
time > the age of the remnant is also indicated by the large size of the X-ray nebula (Seward et 
al. 1984b). The lack of rapid cooling can account for the lower efficiency of X-ray production. 
Gil. 2-0. 3 may also have a relatively flat spectrum (Vasisht et al. 1996), but there is insufficient 
information to draw firm conclusions on this source. 

The observed spectral index of the Crab Nebula is consistent with injection at the shock front 
with p = 2.2. This value of p is consistent with expectations for particle acceleration in a highly 
relativistic shock front; Bednarz & Ostrowski (1998) found that p — > 2.2 in the limit of a high 
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Lorentz factor shock wave. However, N157B does have a steeper spectral index, which suggests 
p ~ 3, and p = 2.2 does not appear to apply to some of the other objects. In the context of the 
present model, detailed spectral study of X-ray pulsar nebulae will be valuable for the determination 
of the particle spectrum created in a relativistic shock front. The same uncertainty is present in 
models of GRB (gamma-ray burst) afterglows. The model presentation here was chosen to parallel 
that for GRB afterglows (e.g., Sari et al. 1998) in order to show the similarities. 

The present models can be used to predict the spin-down power of the putative pulsars in 
wind nebulae where a pulsar has not yet been observed. The value of the models over a purely 
empirical L x — E relation is that additional information about the nebula should allow a more 
precise determination. An excellent example is provided by the Chandra observations of G21.5-0.9 
(Slane et al. 2000). Slane et al. find a photon index ax ~ 1-5 for an inner core and ax ~ 2.0 for 
the whole wind nebula. The spectral steepening is that expected for synchrotron burn-off. A model 
like that for the Crab Nebula should apply, with high efficiency of conversion of spin-down power to 
X-ray emission. For a distance of 5 kpc, the observed X-ray luminosity is 2.1 x 10 35 erg s^ 1 (Slane 
et al. 2000), from which I estimate E « 5 x 10 36 erg s -1 , as compared to E « 3.5 x 10 37 erg s _1 
deduced by Slane et al. (2000). Detection of a pulsar will allow a test of this prediction. 

Support for this work was provided in part by NASA grant NAG5-8130. 
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Table 1. Properties of Pulsars and their X-ray Nebulae 



PSR Supernova Distance E Photon L x L x /E Refs. a 



Remnant (kpc) (ergs s 1 ) Index (ergs s 1 ) 



J0537-69 


N157B 


50 


4.8 


X 


10 38 


2.5 


5.6 x 10 36 


0.01 


(1),(2) 


B0531+21 


Crab 


2 


4.7 


X 


10 38 


2.10 ±0.01 


2.3 x 10 37 


0.05 


(3),(4) 


B0540-69 


N158A 


50 


1.5 


X 


10 38 


2.04±0.06 


8 x 10 36 


0.05 


(2),(4) 


B 1509-58 


MSH 15-52 


4.2 


1.8 


X 


10 37 


1.90 ± 0.04 


1.5 x 10 35 


0.01 


(5),(6) 


B0833-45 


Vela 


0.25 


7.1 


X 


10 36 


1.67 ±0.01 


1.7 x 10 33 


0.0002 


(6),(7) 


J1811-19 


Gil. 2-0.3 


5 


6.4 


X 


10 36 


1 /I +0-54 i on+0.28 
i - 4 -0.68' J - 8U -0.54 


4.9 x 10 34 


0.008 


(8),(9) 


B1951+32 


CTB 80 


2.5 


3.8 


X 


10 36 


1.8 ±0.1 


6 x 10 33 


0.002 


(6),(10) 


B1853+01 


W44 


3 


4.3 


X 


10 35 


2 

z -°-0.9 


4 +30 x 1Q 33 


0.009 


(11) 



"The references to the data are: (1) Marshall et al. 1998; (2) Wang & Gotthelf 1998; (3) Pravdo & Serlemitsos 1981; 
(4) Seward, Harnden, & Helfand 1984a; (5) Seward et al. 1984b; (6) Kawai et al. (1998); (7) Harnden et al. (1985); 
(8) Torii et al 1999; (9) Vasisht et al. 1996; (10) Safi-Harb, Ogelman, & Finley 1995; (11) Harms et al. 1996 
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